7 February 2008 



UNIL-IPT-00-10 



The role of the input in Roy's equations for 

TT-TT scattering 

G. Wanders 

Institut de Physique Theorique, Universite de Lausanne 
CH-1015 Lausanne, Switzerland 
e-mail : ger ar d . wanders @ipt . unil. ch 

Abstract 

The Roy equations determine the S- and P-wave tt-tt phase shifts on 
a low energy interval. They allow the derivation of threshold parameters 
from experimental data. We examine the solutions of these equations 
that are in the neighborhood of a given solution by means of a lineariza- 
tion procedure. An updated survey of known results on the dimension 
of the manifold of solutions is presented. The solution is unique for a 
low energy interval with upper end at 800 MeV. We determine its re- 
sponse to small variations of the input: S-wave scattering lengths and 
absorptive parts above 800 MeV. We confirm the existence of a univer- 
sal curve of solutions in the plane of the S-wave scattering lengths and 
provide a control of the decrease of the influence of the input absorptive 
parts with increasing energy. A general result on the suppression of 
unphysical singularities at the upper end of the low energy interval is 
established and its practical implications are discussed. 
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1 Introduction 



Low-energy tt-tt scattering is a major testing ground of chiral perturbation the- 
ory 0] . Some of its coupling constants are directly related to the tt-tt threshold 
parameters. At present this relation is established at the two-loop level Jl], 
As it is impossible to measure tt-tt scattering at threshold, this relation cannot 
be exploited directly. A reliable extrapolation of the available experimen- 
tal data down to threshold is required. Such an extrapolation is performed 
presently P] with the aid of the Roy equations @, H. These equations are 
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based on the analyticity, the crossing symmetry and the unitarity of the tt-tt 
partial wave amphtudes. The S- and P-wave Roy equations are solved in by 
means of elaborate numerical methods. In this work we discuss aspects of the 
problem which allow an analytical approach and our effort is complementary 
to the work in 0. 

The Roy equations contain as input the S-wave scattering lengths, the S- and 
P-wave absorptive parts above an energy Eq ( which will be called the "match- 
ing point" ) and driving terms coming from the higher partial waves. The Roy 
equations determine, at fixed elasticities, the S- and P-wave phase shifts below 
the matching point. This is a difficult, non-linear, problem that cannot be 
solved analytically. Here we restrict ourselves to questions that can be an- 
swered by linearization and which allow a partly analytic treatment. These 
concern the multiplicity of the solution and its sensitivity to small variations 
of the input. 

Such questions have already been treated in |^ in conjunction with the early 
phenomenological applications of the Roy equations and we can reduce the 
discussion of the multiplicity of the solution to the statement of our old results. 
The matching point used in [0 is at 1.13 GeV whereas the one used nowadays 
in P] is at Eo = 800 MeV. The answers to our questions depend on the choice 
of Eq. The solution is non-unique if Eq = 1.13 GeV and becomes unique when 
Eq = 800 MeV. The response to variations of the input also depends strongly 
on the position of the matching point and our analysis in has to be updated. 

The Roy equations with arbitrarily chosen input make up a well defined math- 
ematical problem. A peculiar feature of this problem is that its solutions 
exhibit unphysical singularities at the matching point [we exclude throughout 
matching points coinciding with an inelastic threshold]. The physical inputQ 
is therefore a special one admitting at least one solution, the physical phase 
shifts, that is regular at Eq. Inputs with solutions regular at Eq have been 
called "analytic inputs" in in the context of simplified elastic one-channel 
Roy equations. The discussion of that class of inputs is extended here to the 
case of the complete coupled inelastic Roy equations. The main conclusion is 
that an analytic input admits a unique solution that is regular at the matching 
point. The non-uniqueness problem is thus evacuated. 

Although non-uniqueness and singularities at Eq are physically excluded, they 
show up in practical calculations because one is working with an approximate 
input which is not exactly an analytic one. An arbitrary variation of an analytic 

^By physical input we denote the input corresponding to scattering amphtudes which 
would be measured in the absence of isospin violation. Ideally, our physical amplitudes are 
those provided by QCD; in practice they are given by the available analyses of experimental 
data assuming isospin symmetry. 
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input produces a non-analytic one and induces singularities at Eq even if the 
choice of Eq guarantees uniqueness. 



We find that one may stay close to an analytic input by correlating suitably 
the variations of two distinct pieces of the input. This comes mainly from the 
matching point at Eq = 800 MeV that is near the p-meson mass. For instance, 
singularities at Eq are largely suppressed by correlating variations of the isospin 
and 2 S-wave scattering lengths Oq and Oq. This confirms the existence of 
a physically acceptable family of solutions along a "universal curve" in the 
(oq, aQ)-plane Similar suppressions of singularities take place if a localized 



variation of an input absorptive part is combined with a variation of one of the 
scattering lengths, Oq for instance. The response to such variations provides 
information on the sensitivity of the phase shifts to the uncertainties on the 
input absorptive parts. We find a very weak sensitivity to the uncertainties 
above 1 GeV. All our results are in qualitative and quantitative agreement 
with those obtained numerically in [Q. 

The coupling between the S- and P-wave channels built into the Roy equations 
is a manifestation of crossing symmetry. The practical implications of this 
symmetry are not well understood and the effects of variations of the input 
might be expected to provide some insight. We find that this is not really the 
case. In our framework the response to a change of the input absorptive part 
in one channel is largest in the same channel but the responses in the other 
channels are not much smaller. All we may say is that crossing symmetry 
produces a substantial coupling of the three S- and P-wave channels, but we 
do not recognize very striking features. 

The paper is organized as follows. The linearization procedure developed 
in 1^, 1^ is described, and the status of the uniqueness problem is outlined, 
in Section 2. Section 3 is devoted to the response to variations of the S- 
wave scattering lengths and the existence of a universal curve. The effects 
of correlated localized variations of input absorptive parts and variations of a 
scattering length are presented in Section 4. Variations of the driving terms 
are also briefly discussed in that section and our conclusions are displayed in 
Section 5. The fact that an analytic input admits only one solution that is 
regular at the matching point is a crucial result. We find it convenient to 
separate its proof from the presentation of phenomenological results and to 
explain it in Appendix A. Our approximation scheme for the determination of 
linear responses is described in Appendix B and Appendix C gives a list of the 
kernels entering the S- and P-wave Roy equations. 
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2 Solution manifold of the S- and P-wave Roy 
equations 



To set the stage we recall the main features of the S- and P-wave Roy equa- 
tions [^]. They relate the real and imaginary parts of the S- and P-wave tt-tt 
scattering amplitudes at low energies, below the matching point Eq : 



Re/,(s) 



■4)E 

j=0 



so 

dx — 

4 X 



5i 



X 



+ Rij{s,x) 



lmf^ix) + Ms), (2.1) 



2 = 0,1,2. To lighten the writing, our notation differs from the standard 
one: /o and /2 are the isospin 1 = and 1 = 2 S-wave amplitudes and /i 
is the isospin 1 = 1 P-wave. We return to the conventional notation // in 
the presentation of final results. The variables s and x are squared total CM 
energies in units of (M^r = pion mass, sq = {Eo/M-^Y). The equations ( p.l| ) 
contain singular diagonal Cauchy kernels and regular kernels Rij which are 
displayed in Appendix C. 



The (pi are input functions 

0j(s) = + (s - 4)|ci(2ao - 602) 



(2.2) 



-E 

j=0 



TT 



dx- 



+ Rij{s,x) 



Aj{x) + ipi{s)y 



In this equation oq and 02 are the isospin and 2 S-wave scattering lengths, 
ai = here and 

_ 1 _ 1 _ 1 _ 
^° ~ 12' ~ 72' ~ ~24' 
the Ai are the absorptive parts above the matching point : 

Ai{s) = lmfi{s), s > So, 



(2.3) 



(2.4) 



and the ipi are so-called driving terms describing the contributions of the higher 
partial waves (/ > 2). They have partial wave expansions converging in [4, sq] 
as long as Sq < 125.31 The equations ( ^TTP constrain the S- and P-waves 
on [4, So] at given input (a^, Aj, ^/'j). Unitarity implies 



Ms) 



2ia(s] 



2iSi{s) 



'1--, 



(2.5) 



where 6i is the channel i phase shift and rji is the elasticity parameter (0 < 
rji < 1) which we incorporate into the input. 



At given input the equations ( |2.1| ) are coupled non-linear integral equations for 
the phase shifts 6i on the interval [4, sq]. To be acceptable, a solution of these 
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equations has to provide absorptive parts below sq that join continuously the 
inputs Ai at that point: 



lim — -— 

s/'so 2a{s) 



[1 - r]i{s) cos(25i(s))) = ^^(so). 



(2.6) 



This boundary condition has to be added to the equations 



The Roy equations being singular, the uniqueness of their solution is by no 
means guaranteed. We sum up the discussion of that point using the technique 
developed in 0, 0. This technique will be our main tool throughout this article. 

We assume we have a set of phase shifts Si satisfying the equations ( ^.1| ) and 
(|2.6| ), the amplitudes fi being given by ( |2.5| ). We ask if these equations have 
other solutions 5- with the same input. If the 6'^ are infinitesimally close to the 
Si the differences {6^ — 6i) obey the linearized coupled equations 



v-^ 1 f^" 1 
cos{2Si{s))hi{s) = z2~ J '^^ 



+Rij{s,x) 



X 



sm{26j{s))hj{s) 



where 



hiis) 



The boundary conditions ( p. 61) imply 



0, 



(2.7) 
(2.8) 

(2.9) 



i.e. 5,'(so) = Si{so). The homogeneous equations ( |2.7| ) with boundary condi- 
tions (|2.9|) may have non-trivial solutions because of the presence of Cauchy 
kernels. The uniqueness or non-uniqueness of the 6i depends on the existence 
of such solutions. 



If the regular kernels Rij are omitted, the equations ( p.7| ) decouple and one 
recovers the one-channel problem discussed in ^ . The existence of non-trivial 
solutions of this problem depends on the value of the phase shift 6i at the 
matching point sq. We assume that 5i(so) > — 7r/2. There is no solution if 
—n/2 < 6i{so) < Tx 12. If 5i(so) > 7r/2, the general solution is 



where 



with 



GAs] 



So 



So - s 



mi 



exp 



4 X — S 



rrii 



-5i(so) 

TT 



(2.10) 
(2.11) 
(2.12) 
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[x] is the greatest integer smaller than x (as in |p, sq is chosen in such a way 
that Si{so) is not an integral multiple of 7r/2). The last factor Pj in the r.h.s. 
of (|2.1CI|) is an arbitrary polynomial of degree rrii — 1. 



The general solution of the complete set of coupled equations (|2.71 ) has a form 
similar to ( [2.101) : 

hi{s) = {s- A)G,{s) [n{s) + H,{s)] (2.13) 
with corrections Hi M. 



The Pi are again arbitrary poynomials of degree rrii — l: mj is given by ( p.l2|) if 



^i(so) > 7r/2; mj = and = if |5i(so)| < vr/2. The functions Hi are regular 
on [4, So] and are solutions of a set of coupled non-singular integral equations: 

1 r.^r..^..:^....^^.H,{x)-H,{s) 



Sm„oHi{s) / dxGi{x)sm{2Si{x)) 

vr J4 



X — s 



, 1 fso 

dxRi,{s,x)G,{x)sm{26,{x))[P,{x)+H,{x)]. (2.14) 



According to definition ( p.ll|) we have 

G,{s)r^iso-sy^ (2.15) 

2 

for s ~ So with 7j = — 5i(so) — rrii. This shows that Gi vanishes at Sq if 

71 

Si{so) > and diverges at that point if 6i{so) < 0. Due to the regularity of 
Hi at So the boundary condition ( |2.9| ) is automatically fulfilled if (5j(so) > 0. 
If — 7r/2 < (5j(so) < 0, Hi has to vanish at Sq (remember that Pj = in this 
case). 

We now apply these results to the uniqueness problem of the physical vr — vr S- 
and P-waves as solutions of the Roy equations (pTl|). The input {ai, Ai,ipi,rii) 
is identified with the physical one and we take the physical phase shifts as our 
master solution 6i of the Roy equations. The physical isospin S-wave and 
isospin 1 P-wave phase shifts being positive |12|, we have mo > 0, mi > 



and 7o > 0, 7i > 0. On the other hand, the isospin 2 S-wave phase shift 62 
is negative (— 7r/2 < ^2 < 0), m2 = 0, P2 = and 72 < 0. Consequently, the 
three boundary conditions (|2.9| ) are satisfied if 



H^iso) = 0. (2.16) 
As solutions of the equations (p^.l4D, the Hi are linear functionals of the poly- 



nomials Pq and Pi. Condition (|2.16| ) gives a homogeneous linear equation 
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relating the coefficients of these polynomials and reduces by one the number 
of free parameters. If mQ + rrii > 1 we are left with mo + mi — 1 free parameters. 
There is no non-trivial solution if mo + rrii < 1. If rriQ + mi > 1, the physical 
phase shifts are embedded in a d-dimensional manifold of solutions of the Roy 
equations with d = mo + mi — 1. If mo + mi < 1, they form an isolated solution 
of these equations. 

The actual values of mo and mi depend on the choice of the matching point so. 
Taking into account the known behaviour of the physical phase shifts ||12| , one 
finds four different situations when Eq = y/soMj, is lowered from 1.15 GeV to 
threshold. 



1. 1 GeV< Eq < 1.15 GeV In that interval, tt < 6o{so) < 37r/2, 7r/2 < 
(^il-So) < ^- This gives mo = 2, mi = 1 and d = 2. The physical S- 
and P-waves are members of a two-parameter family of solutions of the 
Roy equations at fixed physical input and fixed phase shifts at Sq. The 
physical solution can be selected by imposing the physical values of the 
position and width of the p-meson. 

2. 860 MeV< ^o < 1 GeV We now have 7r/2 < 6i{so) < n, i = 0,1 and 
mo = nil = 1, d = 1. The polynomials Pq and Pi reduce to constants 
related by the equations (p.l6|) . The physical amplitudes belong to a 
one-parameter family of solutions. The position of the p-meson can be 
used as a parameter. 

3. 780 MeV< Ep < 860 MeV In this interval mo = 0, mi = 1 and = 
because < 5o(so) < tt/2, n/2 < Si{so) < n. The polynomial Pq 
vanishes and Pi is a constant which is set equal to zero by condition 
( |2.16| ). The physical amplitudes form an isolated solution of the Roy 
equations. Position and shape of the p-resonance are determined by the 
input. 

4. 280 MeV< Ep < 780 MeV Here < 6i{so) < 7t/2, i = 0,1, mp = mi = 
and both Pp and Pi vanish. The physical amplitudes again define an 
isolated solution. 



The above results concern the mathematical problem defined by equations 
(p.l|) , ( p.5|) and ( p.6|) . Due to the behavior (|2.15|) of the Gi at So, the repre- 
sentation ( |2.13| ) implies that if there are solutions 6^ in the neighborhood of 5i 
they are singular at sp and exhibit cusps at that point. These singularities are 
unphysical because the choice of sp is arbitrary. In fact, all solutions of the Roy 
equations with arbitrary input are singular at sp. The physical amplitudes be- 
ing regular at sp, the physical input has to be such that the corresponding Roy 
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equations have at least one solution which is nonsingular at sq and coincides 
with the physical amplitudes. Among all possible inputs the physical input is 
a very special one: it is an an analytic input in the sense of Ref. It has 
been shown there that in simplified one-channel Roy equations with analytic 
input there is only one solution which is regular at Sq. This crucial result is 
extended to the present realistic case in appendix A. 

We see that there is no non-uniqueness problem when working with the exact 
physical input. For instance in case 1 above, one could vary the position 
and width of the p-resonance in the two-parameter family of solutions. The 
singularities at sq would disappear at the physical values of these parameters. 
In practice, however, the physical input is only known approximately and one 
is not really working with an analytic input. Therefore singularities are present 
at So and non-uniqueness cannot be avoided if Eq > 860 MeV. We have to put 
up with these unpleasant features which are merely consequences of a deficient 
knowledge of the physical input. 

From now on we choose the matching point used in the low-energy extrapola- 
tion based on the Roy equations performed in : Eq = 800 MeV (sq = 33). 
Non-uniqueness is avoided but there are unwanted cusps at the matching point. 
It turns out that some of these cusps are in fact a helpful tool. Their suppres- 
sion provides insights into the correlations constraining the scattering lengths 
of an analytic input. This will be illustrated repeatedly in this paper. 



3 Varying the S-wave scattering lengths: uni- 
versal curve 



We come now to our main topic, the linear response to small variations of 
the input. We proceed along the same lines as in the previous section. Start- 
ing from the solution 6i with input {ai, Ai,ijji,rji), we determine the solution 
5- produced by a slightly modified input in linear approximation. To obtain 
quantitative results, we need a model for the 6i which provides an acceptable 
representation of the physical phase shifts. We use the Schenk parametriza- 
tion [g: 



Si 



tan 



-1 



S - Zi 



where 



f^ 



for i 
for i 



0,2, 



(3.1) 



(3.2) 
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The values of the parameters are given in Table 1 and the phase shifts are 
shown in Fig. 1. 



i 




k 


Q 







0.200 


0.245 


-0.0177 


39.3 


1 


0.035 


2.76-10-^ 


-6.9-10"5 


31.1 


2 


-0.041 


-0.0730 


-3.2-10-^ 


-37.3 



Table 1: Values of the coefficients in the parametrization ( |3.lD of the physical S- 
and P-wave phase shifts. 



The elasticities rji are very close to 1 below our matching point and will be set 
equal to 1 in all our numerical results. 

In the present section, we vary only the S-wave scattering lengths: 

ai ^ a[ = ai + 6ai, i = 0,2. (3-3) 

One of the main goals of the low-energy extrapolation of the experimental data 
lies in the determination of these scattering lengths. This cannot be achieved 
directly by solving the Roy equations because the scattering lengths enter into 
the input of these equations. However, they can be predicted in an indirect 
way because the physical input is an analytic one. Consequently, the scattering 
lengths are not independent of the other pieces of the input. If we know the 
physical Ai, ipi and rji we may solve the Roy equations for arbitrary scattering 
lengths ttj. According to Proposition 1 in Appendix A, their physical values 
are obtained by varying these Oj until one arrives at a solution which is regular 
at sq. In practice, when working with an approximation of the physical Ai, 
ipi and r]i, the scattering lengths have to be varied until the corresponding 
solution of the Roy equations can be declared a good approximation of the 
solution of the problem with exact input. This is precisely the procedure used 
in and it is instructive to have an explicit control of the response to the 
variations fB.3|). 



Our task is to determine functions hi defined as in equation (|2.8|) . In the 
linearized scheme the representation (|2.13|) is replaced by 



h^{s) = G,{s) 1^ + - 4) [(-s)™'c,(25ao - BSa^) +p5,,, + Hi{s)]\ . (3.4) 
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E (GeV) 



Figure 1: The tt-tt S- and P-wave phase shifts according to the Ansatz and data 
points obtained from analyses of experiments: 6q and 5| from [14| and 6q from |15]. 
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We recall that mo = m2 = 0, mi = 1, Pq = -P2 = and Pi is a constant which 
we call p. The functions Hi are the solutions of inhomogeneous extensions of 
equations (|2.14|) . They have the form of equations (|B.1|) with 



2 

j=0 

1 r^o 



dx Rij{s,x)Gj{x) sm{26j{x)) x 



6a^ 



(x-4)G,(4) 



+ [(-x)™^ Q(25ao - 55a2) +p5,-,i] ^ . (3.5) 



The solutions Hi depend linearly on p and this constant is fixed in such a way 
that /i2(so) = 0, i.e. 

5a2 



^2(4) 



+ {so - 4) [c2{26ao - 5^02) + H2{so)] = 0. 



(3.6) 



We refer the reader to Ref. [0] for a derivation of equations (|3.4| ) and (tB.ll) . 



To obtain the /ij we have to evaluate the molulating functions Gi on the right- 
hand side of ( pr^ ) and find the solution Hi of ( |B.lD and (|3.5|). The functions 
Gj defined in ( |2.11| ) and obtained from the model ( p.l|) are shown in Fig. 2. 
The exponents 7^ appearing in ( 2.15 ) are 



70 = 0.89, 71 = 0.19, 



72 



-0.20. 



(3.7) 
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The small value of 71 comes from the fact that sq is close to zi, the position 
of the p-resonance. This leads to the spectacular cusp of Gi at sq, seen in 
Fig. 2. As h2 behaves as [sq — s)G2 at sq, it is this product which is relevant 
and shown in Fig. 2. The exponents 70 and 72 + 1 are close to 1 and the cusps 
of Go and (sq — s)G2 are not visible in the figure. 





12 16 20 24 28 32 



Figure 2: The functions Gi defined in ( 2.11 ) appearing as factors in the responses 

So - s 



to variations of the input. 



Go, 



Gi, 



So 



-G2. 



The Hi are slowly varying and ( P^ ) tells us that hi has a sharp cusp at Sq. 
The 6i defined in are models of the physical phase shifts: they are regular 
at So and meant to be produced by an analytic input (oj, Aj, ^/'j, T^j). We see 
that 6[ obtained from ( |2.8| ) has a sharp cusp at sq. This singular behavior 
is a visible signal that the modified input (a^ + Sai, Ai,'ipi,rii) is no longer an 
analytic one. The differences {6'^ — 6i) are linear in Sqq and 6a2'- 

6l{s) - 6,{s) = Gi{s) [U{s)6ao + fi2{s)6a2] , (3.8) 

where fio and fi2 are regular at so. We see that the cusp of 6'i is suppressed if 

/io(so)5ao + /i2(so)5a2 = 0. (3.9) 

There is a direction in the (oq, a2)-plane along which 5[ has no cusp. There 
is still a singularity at so but no infinite slope. This indicates that the input 
(oj + 6ai, Ai,ipi,rii) is close to an analytic input if condition ( p.9| ) is satisfied. 
An analytic input is not an isolated one: it is transformed into new analytic 
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inputs by suitably correlated variations of its ingredients. Our finding and 
the fact that 6'^ and 62 have no visible cusps show that the physical input is 
transformed into nearly analytic inputs by variations of the scattering lengths 



obeying (3^). One can show that variations of the scattering lengths alone 
cannot transform an analytic input exactly into an analytic one. A movement 
along a direction in the (oq, a2)-plane has to be accompanied by modifications 
of the remaining pieces of the input if one wants to keep it exactly analytic. 
Our results show that these modifications are small and we confirm at the local 
level the existence of a one-parameter family of nearly analytic inputs along a 
universal curve, 02 = 02(09) in the (oq, a2)-plane []TT|. 

To go beyond qualitative results, equations (|B.1|) have to be solved and we 
apply the approximation scheme described in Appendix B. In this scheme, the 
Hi have the form 

H,{s) - (^^y^ {Hi,o{s)6ao + Hi,2{s)^a2 + Hr,^{s)p) , (3.10) 

where the Hi^k are second-degree polynomials. Once we have determined Hi^k 
and when p has been fixed by (|3.6|) , we find that the ratio 5a2/5aQ defined 
in ( ^l9|) , for which the cusp in 5[ is suppressed, is equal to 0.197. In fact the 
ratio /i,o(s)//i,2(s) is nearly constant and equal to its value at sq on the whole 
interval [4, sq]. We identify the ratio 0.197 with the slope of the universal 
curve at (09,02): it coincides with the slope found in 0]. There is a strong 
compensation of the two terms in the right-hand side of equation ( p.8|) when 
« = 1 if one moves along the universal curve. This compensation is maximally 
removed in the orthogonal direction where 5a2/5aQ = —5.08. This is illustrated 
in Fig. 3 which displays the relative phase-shift differences (5- — 5j)||/5j and 
((5j' — Si)j_/6i. The differences (5- — 6i)\\ are obtained when the point (00,02) 
moves along the universal curve 

(5oo = 5o|| cos^^ij, 5o2 = 5o|| sin6'||, (3-11) 

with = tan~^ 0.197 = 11°. The differences (5^ — 5j)_L are obtained in response 
to a displacement {6aQ,Sa2) normal to the universal curve, (5oo and ^02 being 
given by ( p.ll|) with (5o|| replaced by Sa± and 6'|| replaced by 9± = 101°. 
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Figure 3: Relative responses (Sf — Sl)/6l to variations of the S-wave scatter- 
ing lengths, (a) displacement (|3.11| ) along the universal curve, 
(b) displacement orthogonal to the universal curve, 6aj_ = 0.05|aQ|. 



To assess the quality of the results displayed in Fig. 3, Table 2 gives the values 
of the Xi defined in (|B.10|) in the parallel and orthogonal directions. The values 
of the xt S'i'e acceptable whereas those for xf are surprisingly small. 
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i 





1 


2 




3.7-10-^ 
9.5-10-3 


1.9 - 10-3 
3.1 - 10-2 


2.9 - 10"^ 
1.8 - 10-2 



Table 2: Accuracy of the approximate values of {5[ — and {5[ — 5i)_\_. The mean 
relative quadratic discrepancies Xi ^-re defined in (B.IO). 



The relative variations of the phase shifts in the direction of the universal curve 
are decreasing functions of s. The S- waves have peaks at threshold, whose 
sizes are dictated by the values of da^ and(5a2. The pattern in the orthogonal 
direction is different and more complicated. The effects of the variation of the 
scattering length spread over the whole interval [4, sq] and {5[ — 5i)^/5i has a 
cusp which cannot be overlooked. 

The overall size of the variations {5[—5i)_[_ is significantly larger than that of the 
corresponding {5[ — 5i)\\. To characterize this fact quantitatively we evaluate 
the mean values of the absolute ratios over the interval [4, Sq — 2] (sq — 2 instead 
of So as upper limit to avoid effects of the cusp in {5[ — 5i)±) 

_ I iK - 

We take 5a± = Sa\\, i.e. we assume that the same distance is covered along and 
perpendicularly to the universal curve, and find 

Po = 15.7, pi = 217, p2 = 18.6. (3.13) 

These large values reflect the sharp definition of the universal curve obtained 
in §. 



(3.12) 



4 Combined variations of input absorptive parts, 
driving terms and scattering length 

This section is mainly devoted to variations 6Ai of the absorptive parts Ai 
above the matching point. As shown in [0 the response is obtained from the 
following Ansatz for the functions hi in ( |2.8|) : 

h,{s) = G,{s){s - 4)[p5,,i + F,{s) + H,{s)], (4.1) 
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where 

^' Ti Jso x"^- X - Ax - s GAx) ^ ' 



The functions Hi on the right-hand side of eq. ( [4.1|) are solutions of the equa- 
tions ([B.l|) with 



Z,(s) = 5:%(s), (4.3) 

i=o 



where 



%(s) = Yij{s) + - AxR,,{s,x) sin(25,(x))G,(x) [p5,,i + F,(x)] (4.4) 
vr J4 



and 

1 r°° 



I roo I 

^iAs) = -j ^x-—^Rij{s,x)5Aj{x). (4.5) 



Equations (^I^)- (|0| ) have been derived in IQ. For simphcity we assume that 
the 5Ai vanish at Sq and the boundary condition ( p.9|) remains unchanged: 
hiiso) = 0. 

Using the analyticity properties of the kernels Rij referred to in Appendix B, 
the integral in the right-hand side of eq. (|4.4| ) can be transformed to give 



foo I dAiu) 
Z^J{s)=pQ,{s)S,,l+ duM,,{s,u) -r^rri, (4.6) 

where Qi is a known polynomial and 

Mij{s,u) = -— (J}dxRij{s,x)Gj{x) f - j -— -. (4.7) 

The contour V encircles the segment [— (s — 4),0]. Formula ([4.7|) affords an 
explicit evaluation of Mjj once the Gj have been approximated by polynomials, 
as in appendix B. The functions Hi are determined by applying the method of 
that Appendix. The condition /i2(so) = fixes j9 as a linear functional of the 
5Ai and the differences 5'i — 6i resulting from ( p. 81) and ( [4.1|) can be written as 

6l{s)-6,is) = J2 duK,,{s,u)6A,{u). (4.8) 



j=o 

The kernel Kij{u) gives the effect on the channel i phase shift of a variation 
of the channel j absorptive part at point u {u> sq). It follows from ( [4. ID that 
Kij{s,u) is proportional to Gi{s) and Kij exhibits, as a function of s, a sharp 
cusp at s = Sq. An arbitrary variation of the input absorptive parts transforms 
an analytic input into a non-analytic one. We correct this partly and stay in 
the vicinity of an analytic input by modifying simultaneously the scattering 
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length a2 and choosing 5a2 in such a way that the cusp in Kij is suppressed. 
This variation 6a2 is a hnear functional of the SAf. 

2 I'oo 

^02 = ^/ duKi{u)5Ai{u) (4.9) 



i=0 ■ 



and Section 3 tells us that equations of the form ( [4. 81 ) are still valid if the Kij 
are replaced by new kernels Kij. 



The kernels Kij and Kij have been evaluated as functions of s for 3 values of u: 
ui = 35 {El = 828 MeV) close to the matching point, U2 = 51 {E2 = 1 GeV) 
and U3 = 100 {E3 = 1.4 GeV). 
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Figure 4: The kernels Kn and Kn at u = 35 as functions of s: Kn includes the 
effect of a variation of Oq suppressing the cusp in Kn. Ku, Kn. 



The passage from Ku to Ku aX u = ui = 35, slightly above the matching point, 
is illustrated in Fig. 4. As must be the case, the large cusp in Ku has disap- 
peared in Ku. The effect of the induced variation of 02 (Ki(ni) = —0.0085) 
dominates the response to the variation of Ai outside the neighborhood of Sq. 
The fact that Ku is larger than Ku is peculiar: Kiq and K12 are much smaller 
than KiQ and K12. 

The values of the kernels Kij &t u = ui determine the responses to small 
variations 5Aj concentrated around that point. If 5Aj is sufficiently small and 
narrow the equations ( [4.8| ) and ( [4.9| ) give 



5[{s) - 5,{s) ~ ^i^,,(s,Mi)A/l„ 5a2 ~ J] «:,K)AA„ (4.10) 
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with 



AAi = / duSAi{u). 



(4.11) 



The relative phase shift differences produced by such variations of the input 
absorptive parts with corresponding variations of the scattering length 02 are 
displayed in Figs. 5, 6 and 7. The AAi have been chosen in such a way that 
the responses are of the order of a few percent. Our linearization should be 
reliable under these circumstances. To describe the situation in physical terms 
we can imagine that the AAj are produced by the insertion of fictitious narrow 
elastic resonances of width Tj at ui. The values of the AAj used in Figs. 5, 
6 and 7 correspond to Fq = 0.76 McV, Fi = 1.56 MeV, F2 = 0.92 MeV. The 
effects of these resonances sitting just above the matching point spread over the 
whole interval [4, sq]. The induced variation of 02 produces a modest peaking 
of (^2 — ^2) / ^2 at threshold. The responses to a variation AAq in the isospin 
S-wave are globally smaller than the effects of variations AAi and AA2 of 
the same size in the other channels. A variation AAi in channel i produces 
a response in the same channel that is enhanced near sq and dominates the 
responses in the other channels. This dominance is significant but not very 
strong in the case of AA2. Apart from these observations we do not discover 
any striking feature characterizeing qualitatively the coupling of the S- and 
P- waves. 



0.003 




-0.012 I ' ' ' ' ' ' ' 1 

4 8 12 16 20 24 28 32 36 



Figure 5: Relative responses ((5/' — Si)/6l to a variation of Aq concentrated on 
u = 35, A^o = -0.1. {5'o'-S'o)/S'o, i^l' -^D/^l {S'o -Si)/6l 
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-0.04 



8 12 16 20 24 28 32 36 



Figure 6: Relative responses {Sf — 6i)/Sf to a variation of A\ concentrated on 
u = 35, A^} = -0.1. (<5§' - <5§)/<5§, {SI' - dl)/dl {Si' - <5g)/(5g. 
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0.03 - 



0.02 - 
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Figure 7: Relative responses (dj' — <^f)/^/ to a variation of Aq concentrated on 

u = 35, AAl = 0.1. (58' - 5§)/<58, ('^i' - ^D/^l (^o'' - 

The Kij are decreasing functions of u without significant change in their shape 
as functions of s. The decrease is rapid just above the matching point. For 
instance, the Kio are scaled down at u — 36 to 70% of their values ai u — 35. 
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To characterize the decrease of the responses when variations AAj are shifted 
to higher energies, we compute averages pij of the absolute values of the relative 
phase shift differences at ui = 35, U2 = 51 and = 100. According to ( [4.10|) 
these are given by 



PijiUk) 



So -4 



ds 



5i 



(4.12) 



Approximate values of the pij{ui) are given in Table 3. The ratios Pij{u2) / Pijiui] 
and Pij{u^) I Pij{ui) show the decrease of the responses at higher energies. None 
of the mean responses to variations located at U2 exceed 11% of the correspond- 
ing responses at Ui. This percentage is reduced to 1.2% when U2 is replaced 
by Ms. Table 4 gives the values of the variations 5a2 coming from ( 4.10|) . 





pijM 


Pij{u2) 
Pij{ui) 


P^i(^3) 
pij{Ul) 


(0,0) 


7.5 ■ 10-3 


0.061 


0.0049 


(1,0) 


4.2 ■ 10-^ 


0.10 


0.010 


(2,0) 


1.1 ■ 10-3 


0.10 


0.010 


(0,1) 


2.4-10-3 


0.10 


0.012 


(1,1) 


2.6 • 10-2 


0.073 


0.0069 


(2,1) 


6.8-10-3 


0.11 


0.012 


(0,2) 


8.3-10-3 


0.080 


0.0069 


(1,2) 


7.5 - 10-3 


0.079 


0.0069 


(2,2) 


2.2 - 10-2 


0.052 


0.0035 



Table 3: Mean relative responses pij{ui) defined in ( [4.12 ) for |A^j| = 0.1 and ratios 



of mean relative responses at U2 and vs. responses at ui, 
Ju^M^ = l GeV, JuiM^ = 1.4 GeV. 



ulM^ = 828 MeV, 



We conclude that the solution of the Roy equations is quite insensitive to the 
errors on the input absorptive parts above = y/usMj, = 1.4 GeV. The 
solution of the Roy equations is most sensitive to the input absorptive parts 
close to the matching point. According to Table 4, the uncertainty in 02, 
associated in our scheme with an error on the absorptive parts at u^, is less 
than 1% of the uncertainty due to the same error at ui. 
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J 


50-2(^1) 


502(^2) 


6a2{us) 




6a2{ui) 


5a2{ui) 





-4.8 ■ 10-3 


0.09 


0.009 


1 


-2.1 ■ 10-2 


0.10 


0.011 


2 


3.6 • 10-2 


0.07 


0.006 



Table 4: Relative variations of the scattering length 02 induced according to ( [4.10 ) 
by variations of the input absorptive parts Aj at ui, \AAj\ = 0.1 and ratios of 
variations at U2 and U3 vs. variations at ui, ^/uiMt^ = 828 MeV, ^/wiMt^ = 1 GeV, 
y/uiM^ = 1.4 GeV. 



We close the discussion of variations of the input absorptive parts with an 
assessment of the accuracy of our resuhs. The errors come from our functions 
Hi. These form an approximate solution of equations ( |B.1D with inhomoge- 
neous terms Zi containing a component ([4.3|) coming from variations of the 



Aj at Uk and a component ( |3.5|) due to the corresponding variation of 02. Let 



Xji^k) be the total discrepancy between left- and right-hand sides of eq. ( |B.1| ) 
defined in (|B.11|) . These quantities are listed in Table 5. All equations (|B.1|) 
are verified at least at the percent level, which is sufficient for our purpose. 





Ui 


U2 


Us 


Xo 


0.012 


0.012 


0.005 


Xi 


0.024 


0.057 


0.016 


X2 


0.017 


0.019 


0.014 



Table 5: Total discrepancies Xj defined in ( B.ll ) to variations of the input absorptive 
part Aj at and the correlated variation of the scattering length 02, ^/ui M-,^ = 
828 MeV, y/H^M^ = 1 GeV, Ju^M^ = 1.4 GeV. 



We close this section with a survey of the response to variations of the driving 
terms ipi in equation (|2.2| ). The tpi are small and approximated by polynomials 
on [4, So] in We consider variations of these polynomials. The Ansatz for 
the functions hi defined in ( |2.8| ) becomes 

h,{s) = {s- A)Gi{s) {p6,,, + H,{s)) , (4.13) 
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where p is a constant and the Hi form a solution of the equations ( [B.l| ) with 

Zi{s) = Sipiis) + Si^iP J dxRii{s,x) sm{25i{x))Gi{x). (4.14) 

As before, the variations of the driving terms are combined with variations 
of 02 such that hi has no cusp at Sq. The resuh shows that large relative 
variations of the driving terms affect only weakly the phase shifts below sq. 
For instance, a reduction of the size of ipo or ip2 by 50% changes the Si by less 
than 5%. In the case of a 50% reduction of ipi the response is smaller than 
0.5%. 



5 Summary and conclusions 

We have developed an approximation scheme to determine the linear response 
of the solution of the S- and P-wave Roy equations with matching point sq = 33 
to small variations of their input (S-wave scattering lengths, S- and P-wave 
absorptive parts above sq, and driving terms). Our results are precise at the 
percent level, which is sufficient for a qualitative insight. Our problem has been 
solved long ago, in a different way, for a higher matching point sq = 70 in [0. 
At So = 33 the solution of the Roy equations is unique, entirely determined by 
their input. 

An arbitrary input leads to a solution that is singular at sq. As the physical 
amplitudes are regular at sq, the physical input belongs to the restricted class 
of our analytic inputs producing a solution that is non-singular at sq. We 
prove that under legitimate assumptions an analytic input has in fact only one 
solution regular at Sq (Appendix A). 

An arbitrary variation of the input transforms an analytic input into a non- 
analytic one and induces responses that are singular at sq- Due to the fact that 
our So is close to {Mp = p-meson mass), the sharpest singularities show 
up as cusps in the isospin 1 P-wave responses. These cusps are suppressed by 
correlating suitably the variations of two pieces of the input. We choose to 
associate in this way variations of the isospin 2 S-wave scattering length Oq to 
arbitrary variations of other components of the input. 

It is instructive to compare our strategy with the procedure used in when 
solving the Roy equations themselves. The solution is parametrized in |^ 
by an Ansatz that is regular at so. As one is working with a non-analytic 
approximation of the physical input, the solution is singular at Sq and cannot 
be fitted exactly by the Ansatz. An approximate solution is constructed by 
a least square procedure tuning simultaneously the parameters in the Ansatz 
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and the scattering length Og in the input. In this way the input is brought close 
to an analytic one and the Ansatz gives a model of the corresponding solution. 
In some of its features this machinery resembles our simple strategy. In fact, 
their equivalence for the computation of responses to small variations of the 
input has been checked in the case of the variation of the isospin S-wave 
absorptive parts displayed in Fig. 5. The response obtained by solving the full 
Roy equations coincides with our result within a few percents. This confirms 
that the main factor tuning Oq in is the avoidance of a cusp in the isospin 1 
P-wave phase shift. 

Our technique shows that one stays in the vicinity of an analytic input when 
moving infinitesimally along a given direction of the (oq, Oq) plane without 
changing the other pieces of the input. This confirms the existence of a so- 
called universal curve at the linear response level. 

We have determined the response to localized variations of the input absorptive 
parts above the matching point. It spreads over the whole interval [4, sq] and 
illustrates the intricate coupling of the S- and P-waves produced by crossing 
symmetry. It shows that the sensitivity to the errors in the input absorptive 
parts decreases rapidly with increasing energy. 

Acknowledgments I thank G. Colangelo, J. Gasser and H. Leutwyler for 
their continual interest in the present work. I am especially grateful to J. Gasser 
for his help in the preparation of the figures. 



Appendix A. Analytic input and uniqueness 

An analytic input {ai, Ai,ilJi,7]i) is defined as an input admitting at least one 
solution of the Roy equations which is regular at the matching point. A precise 
definition is given below. In any case it is an indirect definition: as shown at 
the end of this appendix, we know how to construct analytic inputs but we 
are unable to identify an analytic input by direct inspection. Its components 
are correlated: in particular, the scattering lengths depend on the Ai, ipi and 
rji. Analytic inputs are relevant objects because the physical input belongs to 
that class. The aim of this appendix is to prove that an analytic input has 
only one solution regular at Sq. The requirement of regularity at Sq eliminates 
in principle the uniqueness problem. This result has already been established 
in 01 for simplified one-channel Roy equations. 

To establish our result we need general analyticity properties of the partial 
wave amplitudes /j. Let / be one of them. It is known to be the boundary 
value of an analytic function F on the interval [4, 125.31] This function is 
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holomorphic in a complex domain A extending on the real axis from sl = —28 
to sr = 125.31 and provided with a left-hand cut [sl,0] and a right-hand cut 
[4, sr]. We have 

f{s)=\imF{s + ie), se[4,SR]. (A.l) 



Our matching point Sq being above the first inelastic threshold ii = 16, we 
need properties characterizing the elasticity parameters rj which enter into an 
analytic input. According to (|2.5|), r] is equal to the modulus of the S-matrix 
element (1 + 2io"/) which is the boundary value of 




S{z) = l-2\I^Fiz). (A.2) 

This function is regular in A. Using the relation S{z) = S{z) we write 

r]'^{s) = lim S{s + ie)S{s + ie) = lim S{s + ie)S{s — ie), s G [ii, sr]. (A. 3) 

Although it cannot be derived from first principles [|l^], it is legitimate to 
assume that the inelastic thresholds ik {k = 1,2, . . .) are the only singularities 
of / on [4, sr] and that S has an analytic continuation Sjj into the sheet 
reached by crossing the cut [4, sr] from below between two successive inelastic 
thresholds ik and ik+i [Su depends on the pair {ik, ik+i)]- Equation (|A.3|) gives 

r]'^{s) = \imS{s + ie)Sii{s + ie) s E {ik,ik+i), (A.4) 
as long as ik+i < sr. 

We assume that Sjj is regular in the upper half-plane, in a neighborhood D 
of the segment {ik, ik+i)- Equation (|A.4|) tells us that the real- valued rj'^ is the 



boundary value on {ik,ik+i) of a function holomorphic in D and we apply the 
following general result. 



Lemma 1 Let w be a real-valued function defined on the interval {ik,ik+i)- If 
w is the boundary value of an analytic funtion W holomorphic in D, it is the 
restriction to {ik,ik+i) of a function regular in the domain DU D where D is 
the mirror domain of D: D = {z\z G D}. 



A proof of this Lemma is given at the end of this appendix. It implies that r]'^ 
has an analytic continuation regular in a complex neighborhood of {ik,ik+i)- 
We assume that the possible complex zeros of r/^ are at a finite distance from 
{ik,ik+i)- We choose D sufficiently narrow so that r]"^ is non- vanishing on D 
and we have 
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Lemma 2 // the above conditions are fulfilled rj has a holomorphic continua- 
tion from each interval {ik,ik+i) with z^+i < sr into a complex neighborhood 
of that interval with i^ and ik+i on its boundary. 



We turn now to properties of the full amplitude / and establish 



Lemma 3 The real and imaginary parts of f are separately holomorphic in a 
complex neighborhood of each interval {ik,ik+i) (ik+i < sb,) with ik and ik+i 
on its boundary. Here = 0, 1, 2, . . . with iq = 4. 



This is a well known result in the case of the interval [4, zij For any 

interval we define the function 

ia 1 + r/ + 2ia/ ^ ^ 

on {ik, ik+i) [77 = 1 on (4, ii)]. According to Lemma 2, V has a regular analytic 
continuation into a domain N in the upper half plane - [ik,ik+i] belongs to 
the boundary of - except for poles at the possible zeros of the denominator. 
Using unitarity 

Im/ = a|/p + ^(l-r^2) (A.6) 

we find that Im y = on {ik,ik+i)- Lemma 1 is easily extended to the 
case of meromorphic functions and one concludes that V has a meromorphic 
continuation into N U N. The definition (|A.5|) gives 



We assume again that the zeros of the denominators are at a finite distance 
from the real axis and discover that Re / and Im / are indeed separately 
holomorphic in a neighborhood of {ik,ik+i) contained in A^ U iV. The phase 
shift 6 is also regular in a neighborhood of each {ik,ik+i)- 

We close our preliminaries with the structure of / at an inelastic threshold ik 
with square root singularity. There are four functions a, b, c and d that are 
regular in a circle Ck with center ik and radius p such that 



?7>(sj = exp 



-2 (a(s) + ^/s-ikb{s)) 



6y{s) = c{s) + y/s - ik d{s) 



for s G {ik,ik+i) '■ V> ^md 5> designate respectively the elasticity parameter 
and the phase shift above ik- The amplitude / can be written on {ik,ik+i) in 
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terms of a complex phase shift 5> as 



/ 



c(s) + ia{s) + a/s - ik {d{s) + 



(A.9) 



The value of / below ik, on {ik—p, ik), is obtained through analytic continuation 
of the expression ( [A .91) in the upper half-plane along curves contained in Ck- 
The outcome is determined by a complex phase shift 



S^{s) = c{s) + ia{s) - Vii^ {h{s) - id{s)). (A.IO) 
The elasticity parameter rj^ and the phase shift 5< below ik are given by 



exp 



-2 (^a(s) + \/ik — s d{i 



(5<(s) = c(s) - V4 - 



(A.ll^ 



on [ik — p, ik)- The functions b and d interchange their roles when we cross ik- 
Below ii, 7]^ is equal to 1. This implies 



d = 



(A.12) 



in the case k = 1. 

We summarize our findings in 



Lemma 4 The structure of f at a square root inelastic threshold ik is described 
by formulas (A. 8) and l \A.ll\ ). Equation (A.l^) holds atii. 



After this lengthy preparation we are ready for a complete definition of an 
analytic input. 

Definition 1 An analytic input (oj, Aj, V'i, ^71) contains elasticities fulfilling 
Lemma 2 on (4, sq). It admits at least one solution /j, i = 0,1,2, of the S- 
and P-wave Roy equations with Re /j regular at Sq in the sense that they are 
holomorphic in a circle Cgg : |,s — so| < e. These fi satisfy Lemmas 3 and 4. 

We establish the following 



Proposition 1 Let fi, z = 0, 1, 2, form a solution of the S- and P-wave Roy 
equations with analytic input {ai, Ai,ipi,rii) that is regular at sq and verifies 
Lemmas 3 and 4- A second solution f- of these equations, // 7^ fi, is singular 
at So- 
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This proposition is an extension of Proposition 4 in to the reahstic situation. 
In the following proof we assume ii < Sq < i2, which is true for our sq = 33. 



The proof of Proposition |I| is based on Lemmas 2, 3 and 4. To make sure that 
the Roy equations (|2.1| ) guarantee the required analyticity of the we rewrite 
these equations as follows 

I poo I 

Re <l>i(s) = (s-4)-i dx- — rlm$,(x), (A.13) 

TT J 4 [X — A){X — S) 



where 

2 



^j(s) = Ms)-ai-{s-4:) I Q(2ao - 5a2) + X! ~ dx Rij{s,x)lm fj{x) + ipi{s) | . 

(A.14) 



j=0 



In these equations Im fi{s) = Ai{s) for s > sq, and the ipi are the driving terms 
appearing in ( |2.2| ). The fact that Im $j = Im fi on [4, oo) has been used. The 
equations (|A. 13|) ensure that the $j are boundary values of analytic functions 
holomorphic in C\[4, oo). For x G [4, oo) the kernels Rij{s, x) are holomorphic 
functions of s in C\(— oo,0] and the driving terms are regular in the domain 
A without right-hand cut UTOf. Taking all this into account, ( A.14 ) provides 



a representation of the fi ensuring that they are indeed boundary values of 
functions, holomorphic in the domain A, with right- and left-hand cuts. The 
same conclusion holds for the 

To establish Proposition 1, we show that the // have to coincide with the fi if 
the Re // are regular at Sq. Inversion of the dispersion relations (|A.13| ) gives 



Im $,(s) = Im f,{s) = -(s - 4)i / Re$,(a;) ^ ^ ^ 

VTJIRX — 4 X — S 

The regularity of Re fi at sq implies that Re is regular at that point and 
it follows from ( |A.15| ) that Im is holomorphic in Csg- The relation ( [A. 15] ) 



holds true for $• defined by // and the assumed regularity of Re // at Sq 
implies the analyticity of Im f- in C^q. According to Definition 1, Im fi and 
Im // have analytic continuations holomorphic in a complex neighborhood A^' 
of the interval [s^sq] shown in Fig. 8. We conclude that Im fi and Im fl are 
holomorphic in A^' U Csq ■ As 

Im /;(s) = Im /,(s) = Ai{s) for s E [sq, Sq + e), (A.16) 

Im // and Im / coincide in A^' U Cg^ and 

Im f,{s) = Im /;(s) for s G (zi, sq]. (A.17) 
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Figure 8: Domains of the complex s-plane used in the proof of Proposition 1. 



To complete our proof, we have to extend the equahty ( |A.16| ) below the in- 
elastic threshold ii. The discussion of Lemma 4 shows that phase shifts are 
required to go through ii. The equality of imaginary parts above ii at fixed r]i 
implies 

6i{s) = ±6i{s) modvr, se{ii,so]. (A.18) 

We show that the minus sign must be rejected. If, for a given i, 6'^ = —6i mod tt, 
(|A.8|) gives 

S'i>{s) = c'{s) + ./T^d'is) (A.19) 

with c'{s) = —c{s) modvr, d'{s) = —d{s). According to ( A. 11 ) this would 
produce an elasticity rj'-^ below ii that would differ from the input elasticity 

Vi<- 

In terms of the functions b and c appearing at ii, we now have 

c[ = Ci modvr (A.20) 

whereas 6- = bi is given, rji being a member of the input. Equations (|A.10|) 
and i \KA^ ) give 



^i<i^'^) = - V4 - sbi{s)j modvr, se{ii-p,ii). (A.21) 

Therefore we have 

Im fiis) = Im Ms), s G (^1 - p, Zi). (A.22) 

As both sides of this equations have analytic extensions regular in a neighbor- 
hood N" of (4, ^i) the equality (|A.22|) extends to (4, ii). Thus we have estab- 
lished that and fi have the same imaginary parts on [4, sq] if ii < Sq < 12 
and the Roy equations imply the full equality of these two amplitudes. This 
result extends to arbitrary choices of the matching point. The proof becomes 
easier if sq < ii; it requires more steps if Sq> 12- 



27 



For completeness we prove Lemma 1. We define a function W by W{z) = 
W{z). It is liolomorpliic in tlie mirror domain D and, w being real, we have 
w{s) = lim W{s — ie), s G [ik, ik+i]- We write forzED 

e\0 

W(.) = ±46.^ + ±i6.^. (A.23) 

ZlTT JdD X — Z 2m JdD X — Z 

The first term is the Cauchy representation of W and the second integral 
vanishes because z ^ D. The contributions of the segment [ik^ik+i] to both 
integrals cancel and one is left with the Cauchy representation of a function 
holomorphic in D U D. 

Three remarks close this appendix. 

1. If y^Mjr = 800 MeV, we know, according to Section 2, that the phys- 
ical solution of the Roy equations is an isolated one. The relevance of 
Proposition 1 comes from the possible existence of other solutions with 
^ii^o) = Si{so) + njvr resulting from CDD-pole ambiguities These 
solutions are singular at sq. 

2. The proof of Proposition 1 tells us that the absorptive parts Ai of an 
analytic input are regular on some interval [sq, Sq) above the matching 
point (sg > So + e) and are the analytic continuation of Im /j below sq 
on that interval. The Roy equations ( f^.lD define real parts Re fi above 
sq. On [so, Sq) they are the analytic continuations of the Re fi below 
So- As the interval [sq, Sq) cannot contain an inelastic threshold, all the 
ingredients of the unitarity condition ( |A.6| ) have analytic continuations 



from below sq onto [so,Sq). This implies that (|A.6|) holds on [so,Sq): 
Re fi and Ai are the real and imaginary parts of amplitudes verifying 
unitarity on that interval. This means that they fulfill, at least on [sq, Sq), 
a consistency condition discussed in 0. 

3. Although we have no direct way of checking whether a given input is an 
analytic one, we have a recipe for the construction of such inputs. Take 
a matching point Sq above sq (sq < 125.31) and choose arbitrarily an 
input (a^, A[, ipl, 77-). Let // be a solution of the Roy equations with that 
input, verifying Lemmas 2, 3 and 4. These are expected to be singular 
at Sq but they are regular at sq. Define a new Ansatz (oj, Aj, ■j/'j, r^j) with 
matching point Sq: 

di CI^, 

= i^'i, Vi = v'i on [4, So], 
A'-{s) for s > s'q. 
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The // define a solution /, of tliis new problem, 

Ms) = flis) for 4 < s < So. (A.25) 

This solution is regular at sq and the Ansatz (oj, Ai, ipi, rji) is an analytic 
one. 

Our recipe is of no practical use because it requires the explicit resolu- 
tion of the Roy equations with matching point Sq. The important point 
is that we recognize that an analytic input with matching point sq is 
unconstrained above some Sq, s'q > sq. It is the behavior of the Ai on 
[so, s'q] which is constrained and s'q can be close to sq. 

In our definition an input is analytic with respect to its matching point sq. 
The physical input is special because it generates inputs with matching 
points s'q > Sq {s'q < 125.31) that are analytic with respect to s'q. 



Appendix B. Approximation scheme 

We write the equations we have to solve in Sections 3 and 4 in the following 
way. 

J2X,,[H,]is) = Z,is). (B.l) 
i=o 

Each Xij is a linear and homogeneous functional of the unknown Hj , 

1 r^.n^.^..^ro.f.^^H,{x)-H,{s)\ 



1 j-SO 



X — s 



1 /"^o 

-/ dxRiAs,x)GAx)fiiTi{25Ax))HAx). (B.2) 

TT J4 



The Zi are known functions determined by the variation of the input under 
consideration. The unknown are regular and slowly varying on [4, sq\ and 
we approximate them by polynomials 

N 

^.(s) = 5'"'Ec.,ns". (B.3) 

n=0 

We have to determine the coefficients Cj_„. The X^j become 

N 

where the ^i"^ are known functions obtained by replacing Hj{x) by x^"^^^^^ in 
the right-hand side of Eq. 
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To evaluate these functions we define auxiliary analytic functions Gj, holomor- 
phic in C\[4, sq] 



Gi{z) = ( ) exp 

.So - z. 



2 



TT J4 X — Z 



(B.5) 



They are related to the Gi defined in ( |2.11| ) by their discontinuity Disc Gi 

across the cut [4, sq] , 

^DiscG'i(s) = Gi(s)sin(25i(s)), A < s < sq. (B.6) 

The contribution to X^'^^ coming from the first integral in the right-hand side 
of ( [B.2| ) is transformed into a sum of integrals along a closed contour F sur- 
rounding the segment [4, Sq]: 

-| n— 1 „ n — 1 

- E G^iz)z"^s'^-"^-' = -J2 g-.m-.is''-'^-' (B.7) 

m=0 m=0 

where the (^j^p are the coefficients of the Laurent series of Gi, 

oo -| 

G,{z) = Y.g,,,- (B.8) 

The second integral in the right-hand side of ( p.2| ) is evaluated in a similar 
way by exploiting the analyticity properties of the kernels Rij. At fixed real 
s, s > 4, these are analytic functions of x, holomorphic in C\[— (s — 4),0]. 
Deforming the contour F, we get 

so 

dx Rij{s, x)Gj{x) sin(25j(a;))x" = 

1 .0 _ 

— —7- / dxDisc Rij {s, x) Gj{x)x'^ + polynomial. (B.9) 

2m j-(<i-4) 

The polynomial is determined by the asymptotic behavior in z of the product 
Rij{s, z)Gj{z)z"-. The discontinuity Disc Rij of Rij across [— (s — 4),0] being 
known, we need the Gj on that interval. These smooth functions are approxi- 
mated by third degree polynomials at a level smaller than 1%. This allows the 
explicit evaluation of the integral in the right-hand side of Eq. (B^) and the 
result is a polynomial in s. The xj;^^ are thus approximated by polynomial 
xfj"-* of degree < 6. 

Evaluated along the same lines, the inhomogeneous terms Zi in Eqs. (|B.1|) 
become known functions Zi. According to the equations (p.l| ) the Xi and 
Zi have to be made approximately equal on [4, sq] by adjusting the 3(A^ + 1) 
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coefficients c. 
(p.4| ) and inserting the result into (p.l| )] 



m 



in Eq. ( [B.4| ) [Xi is obtained by substituting X\j for Xl^ 

We keep our calculations simple by 
using polynomials of low degree for the Hi in ( [B.3| ) and choose N = 2. To 
determine the 9 coefficients Q^n, the Xj and Zi are approximated on [4, sq] by 
second degree polynomials using a technique, and these polynomials are set 
equal. This gives 9 equations for the 9 unknowns (in Xj, each xfj^-* is replaced 
by a polynomial of degree 2). The whole procedure is legitimate because the 
Xi and Zi are slowly varying. 



The Xlf and Zi are close to the X^j and Zi, the differences coming only from 
the replacement of the Gi by third degree polynomials on [— (s — 4), 0]. Thus, 
in view of equations ( |B.1[ ), the Cj^„ we obtain must be such that Xi and Zi 
are close to each other on [4, sq]- This can be checked by evaluating the mean 
relative quadratic discrepancies of X^ and Zi 



'X.i 



(so -4)74 



So (Xi{s) 
ds— 



X 



Ziis] 



We can also define a total discrepancy 

"1 2 



Us) 



(B.IO) 



(B.ll) 



The various values we obtain for these quantities are quoted in Sections 3 and 
4. 



Appendix C. The kernels Rij 



Our technique makes extensive use of the analyticity properties of the regular 
kernels Rij in equations ( |2.1| ) and ( ^^.21 ). It is therefore convenient to display 
them explicitly. They are obtained from 4 functions L^, k = 1, . . . ,4: 



Li{s,x) 
L2{s,x) 



1 



x{s - 4) 
1 

x{s - 4) 
1 

x{s-4f { 6 



s-x + 2 + (x- 4) In 1 + 

s - A ' 



s-4 



- X + 2 + (2s + X - 4) In 1 + 

2 ^s - 4 * 



s-4 



- 8s + 4(3x2 - 12x + 4) 



+(2s + X - 4)(x - 4) In 1 + 

s — 4 



s-4 



(C.l) 
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x(s-4)2 i 6 



+ 8s{3x - 1) + 4(3x2 - 12x + 4) 



+(2x + s- 4)(2s + X - 4)(x - 4)^— In (l + - — - 

s — 4 V X 



The Rij are given by 



-Roo('S, 


x) 




x) - 


1 

X 


^02(5, 


x) — 


y ^1(5,2;), 




x) 


- l^s 


x), 




-^22(5, 


x)^ 


^Li{s,x) - 


Roi{s, 


x) 


= 6^2(5, 


x), 




^21 (S, 


x) = 


-3L2{s,x), 


Rio{s, 


x) 


= Ils{s 


x), 




Rl2{s, 


x)^ 


~L3{s,x), 


Rn{s, 


x) 


= 3^4(5, 


x) — 


1 









(C.2) 
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